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Nonequilibrium physics encompasses a broad range of natural and synthetic small-scale systems. 
Optimizing transitions of such systems will be crucial for the development of nanoscale technologies 
and may reveal the physical principles underlying biological processes at the molecular level. Recent 
work has demonstrated that when a thermodynamic system is driven away from equilibrium then 
the space of controllable parameters has a Riemannian geometry induced by a generalized inverse 
diffusion tensor. We derive a simple, compact expression for the inverse diffusion tensor that depends 
solely on equilibrium information for a broad class of potentials. We use this formula to compute 
the minimal dissipation for two model systems relevant to small-scale information processing and 
biological molecular motors. In the first model, we optimally erase a single classical bit of information 
modelled by an overdamped particle in a smooth double-well potential. In the second model, we 
hnd the minimal dissipation of a simple molecular motor model coupled to an optical trap. In both 
models, we find that the minimal dissipation for the optimal protocol of duration r is proportional 
to 1/r, as expected, though the dissipation for the erasure model takes a different form than what 
we found previously for a similar system. 

PACS numbers: 05.70.Ln, 02.40.-k,05.40.-a 


I. INTRODUCTION 

A complete set of principles describing physical phe¬ 
nomena far from equilibrium remains elusive. However, 
considerable progress has been made in the study of 
nonequilibrium processes in recent times. Fluctuation 
theorems relating the probability of an increase to that 
of a comparable decrease in entropy during a finite time 
period have been derived m and experimentally veri¬ 
fied [5H5]. Considerations of Maxwell’s demon and Lan- 
dauer’s principle have led to a better understanding of 
the thermodynamic role of information |1()H12) and new 
fundamental relationships valid for systems far from equi¬ 
librium such as the Jarzynski equality have provided 
deep insights into thermodynamic quantities such as en¬ 
tropy [T!^UT()] . More recently it has been appreciated that 
these nonequilibrium relationships are closely related to 
pioneering work by Bochkov and Kuzovlev HIHIS]. 

Recent work has also shed light on the general prob¬ 
lem of computing optimal protocols that minimize dis¬ 
sipation while driving small-scale systems between sta¬ 
tionary states Optimization schemes for finite¬ 

time thermodynamic processes will be needed for tech¬ 
nological applications in which energetic efficiency is 
paramount |34l I35j . This will be particularly relevant in 
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the decades to come as computational demands approach 
physical limits. 

It may also be the case that evolution sculpted molec¬ 
ular machines such as kinesin and Fq — Fi ATP-ase 
to operate far from equilibrium while maximizing efh- 
ciency [36]. Thus, optimization schemes for nonequilib¬ 
rium transitions also have the potential to unlock phys¬ 
ical principles underlying the function of biological sys¬ 
tems on the molecular level. 

A promising approach utilizes a linear-response frame¬ 
work in which a generalized inverse diffusion tensor in¬ 
duces a Riemannian manifold structure on the space of 
parameters |2B]. Optimal protocols are geodesics of this 
geometry. This idea is developed further in m where 
the machinery of Riemannian geometry is exploited to 
find explicit optimal protocols for a paradigmatic col¬ 
loidal particle model and in m where optimal finite-time 
erasure protocols are computed for a simple classical bit 
model. An extension of this framework to transitions be¬ 
tween nonequilibrium steady states is established in [^ . 

This framework builds on earlier investigations into 
thermodynamic metrics by including the dynamics of the 
driven system. For macroscopic systems, the properties 
of optimal driving processes have been investigated using 
thermodynamic length, a natural measure of the distance 
between pairs of equilibrium thermodynamics states [ 551 - 
I43j . with extensions to microscopic systems involving a 
metric of Fisher information |44l|45|. Slow transitions be¬ 
tween nonequilibrium steady states have also been stud¬ 
ied in terms of thermodynamic metric structure |46j . 
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In [57], the method to calculate the inverse diffusion 
tensor components relied on the (continuous) potential 
being harmonic. In this paper, we demonstrate how the 
inverse diffusion tensor may be computed for a more gen¬ 
eral class of potentials in terms of the equilibrium proba¬ 
bility distribution. Our starting point will be the Fokker- 
Planck equation m, which we assume fully describes 
the physics of the system. Furthermore, we show that 
the inverse diffusion tensor arises naturally through an 
expansion in temporal derivatives in this general setting. 

We use this construction to compute minimal dissipa¬ 
tion for two model systems of physical interest. First, 
we consider a one-dimensional system modeling the stor¬ 
age and erasure of a single classical bit of informa¬ 
tion i37i im 1^ . The erasure of information results 
in energy dissipation according to the Landauer prin¬ 
ciple. Minimizing this dissipation (equivalently, maxi¬ 
mizing erasure efficiency) will likely prove critical to the 
development of future small-scale information process¬ 
ing devices. The first model system consists of an over¬ 
damped colloidal particle diffusing under the influence of 
a continuous double-well potential with a large central 
barrier stabilizing the memory. If the particle is found to 
the left (right) of the origin, the memory value is 1 (0). 
We seek the most efficient protocol altering the shape of 
the confining potential so that the particle will be found 
to the right of the origin with overwhelming probability, 
thus setting the memory value to 0 and erasing the sin¬ 
gle bit of classical information originally encoded by the 
system. 

The erasure cycle consists of a continuous stage in 
which the wells merge and the central barrier is lowered. 
A reset stage in which the potential returns instanta¬ 
neously to its original state and leaves the final probabil¬ 
ity distribution undisturbed completes the erasure cycle 
(see Fig.[^. The inverse diffusion tensor predicts optimal 
erasure cycles in the long duration limit. 

The second model system consists of an overdamped 
colloidal particle diffusing in one-dimension while cou¬ 
pled to a ratchet potential and an optical trap. In this 
simplihed model, the coordinate of the diffusing particle 
may be identified with a mechanical state variable of a 
molecular motor m and the thermal bath consists of the 
huge number of irrelevant degrees of freedom of the liquid 
surrounding the motor as well as the internal degrees of 
freedom of the motor itself and the structures with which 
it interacts. We suppose that our simplistic molecular 
motor couples to an optical trap. The inverse diffusion 
tensor framework predicts the optimal time course for the 
optical trap center, which represents an external drive for 
this simple molecular machine. 


II. THE INVERSE DIFFUSION TENSOR 


For a physical system at equilibrium in contact with 
a thermal bath, the probability distribution over mi¬ 


crostates X is given by the canonical ensemble 

Peq(x, A) = exp/3 [F(A) - E{x, A)] , (1) 

where /3 = (/cbT)”^ is the inverse temperature in natural 
units, F{X) is the free energy, and E{x, A) is the system 
energy as a function of the microstate x and a collection 
of experimentally controllable parameters A. 

In equilibrium, the thermodynamic state of the system 
(the probability distribution over microstates) is com¬ 
pletely specihed by values of the control parameters, but 
out of equilibrium the system’s probability distribution 
over microstates fundamentally depends on the history of 
the control parameters A, which we denote by the con¬ 
trol parameter protocol A. We assume the protocol to 
be sufficiently smooth to be twice-differentiable. 

The average excess power exerted by the external agent 
on the system, over and above the average power on a 
system at equilibrium, is m 


/3(to)'Pex(to) = 


dX^ 

dt 


io 


■ {sx}^ . 


( 2 ) 


Here X = -are the forces conjugate to the control 

parameters A, and <5X(to) = X{tQ) — is the 

deviation of X(to) from its current equilibrium value. 
Applying linear response theory m, 


{SX{to))j 


rio 


Xito - t') ■ [A(to) - A(t')] dt' (3) 


where Xij{i) = ~d^fj*°\t) / dt represents the response of 
conjugate force Xi at time t to a perturbation in control 
parameter at time zero, and 

Ef;‘“)(t) = (5X,(0)5A,(t))^(^^). (4) 

For protocols that vary sufficiently slowly [51] , the result¬ 
ing mean excess power is 


P{to)'Pex{to) 


dy 

dt 


C(A(to)) 


J to 


dX 

dt 


■J to 


for inverse diffusion tensor 

nOO 

Cy= / dt'{6Xj{0)SXy)) 

Jo 

or, more conveniently, 

nOO 

Cij = / dt'{dxj(t){0) dxi(j){t')) 
Jo 


^{to) 


X{to 


( 5 ) 


( 6 ) 

( 7 ) 


where (p{x, A) = — In Peq{x, A). 

We assume our system consists of an overdamped col¬ 
loidal particle obeying the stochastic equation of motion 

X =--d^U{x{t),t) + E{t) (8) 
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for Gaussian white noise F{t) satisfying 


so that 


{Fit))=0, {F{t)F{t')) = ^6{t-t'). (9) 

Here, 7 is the Cartesian friction coefficient and U{x,t) is 
a generic potential growing unbounded as |a::| —>■ 00 at a 
rate to be specified shortly. 

The components in Eq. 0 may be readily computed 
if Eq. Q is linear or, equivalently, the potential is har¬ 
monic |27j . The linearity of Eq. ([^ allows us to write 
general solutions as a linear combination of a homoge¬ 
neous piece dependent only on the initial conditions and 
a particular piece dependent on the Gaussian noise. Such 
a decomposition allows for a straightforward calculation 
of the time correlation functions appearing in Eq. 0 as 
demonstrated in m- 

For more general potentials, such a decomposition is 
impossible and so a different approach must be found to 
compute Eq. 0. We will find the equivalent statistical 
description in terms of the Fokker-Planck equation 

dtp = D[d, m\x,t)p) + dip] = -d^G ( 10 ) 

convenient, where p{x,t) is the position probability den¬ 
sity, G(a;, t) is the probability current, and D is the dif¬ 
fusion coefficient. 

We set out to compute a general expression for the 
inverse diffusion tensor of a driven system obeying over¬ 
damped dynamics. Using the method of Laplace trans¬ 
form |52) . we succeeded in writing the components of 
the tensor entirely in terms of the equilibrium probabil¬ 
ity distribution p^q and cumulative distribution function 
He,: 


C.,(A) 



l9A-neg(3:, A) dxi^eqjx.X) 
Peq{x, A) 


( 11 ) 


We assume that the potential satisfies U{x,\) 00 

as |a;| —>■ 00 . Note that the construction also applies for 
reflecting-wall boundary conditions. 

The inverse diffusion tensor components are given by 



/ OO pOO 

dxo peq{xo,X) dxj4>{xo,\) / dt'niiit'). 
-00 J 0 

(15) 

Note that we have suppressed the dependence of rrii on 
xq and A for convenience. We evaluate dt' ini{t') by 
computing the Laplace transform 

poo 

m,i{s) = / dt' mi(t') (16) 

do 

and taking the limit as s —>■ O'*". 

Integrating by parts, 

dt' ^{t') = s m,(s) - m,(0). (17) 

dt' 

Note that mi(po) vanishes since limt_>oo p(a;, t; a;o) = 
Peq{x,X)', i.e. the system comes to a stationary state 
after a sufficiently long time has elapsed. By definition 
of rrii, 




dx dtp{x,t;xo) dxi(l){x,X). 


In terms of the probability current G{x,t), 

^ TOi(O) - dx dxG{x, s) dxi(j){x, X) 

mi[s) = -— -—--— 


(18) 


(19) 


Therefore, to compute rhi{s), we need the Laplace trans¬ 
form of the probability current. 

The Fokker-Planck equation may be used to derive an 
equation for the probability current: 


dtG{x, t) = D [pU'{x, X)d^G{x, t) + dlG{x, t)]. (20) 


Taking the Laplace transform of both sides, we have 

s G(a:, s) - G(a:, 0) = D[(3U'{x,X) d^G{x,s) + dlG{x,s)\ 

( 21 ) 

which follows from limi_>oo G{x, t) = 0. Multiplying both 
sides by s and defining H{x, s) = s G{x, s), 


where (j){x,X) = — hipeq{x^X) and Peq{x,X) is the equi¬ 
librium distribution. 

We rewrite Eq. as 



' —00 

nOO 


dxo Peq{xo,X) dxi(l)ixo,X)x 

dx p{x,t']xo) dxi(l>{x,X) 


(13) 


where p{x,t;xo) satisfies Eq. ( |Io| ) with initial condition 
p{x,t = 0 ;a;o) = Sx,xo and p{x,t;xo) —>■ 0 for |a;| —>■ 00 . 
For simplicity, define 


s H{x, s)—s G{x, 0) = D[pU'{x, A) dxH{x, s)+dlH{x, s)]. 

( 22 ) 


We may obtain a solution to Eq. (22) by expanding 
H{x, s) as a series in s. If we define lT(x, s) = Hq{x) + 
s F[i(x) -|- ..., then 


Q = pU'{x,X) H'(x)-f H"(x) 


(23) 


Ho{x) - G{x, 0) = D[I3U'{x, A) dxHi{x) + dlHi{x)] 


(24) 


iLfe_i(x) =D[l3U'{x,X) dxHkix) + dlHkix)] (25) 


/ OO 

dx p{x,t-,Xo) dxi4>{x,X) 

-OO 


(14) 


follow from substituting the expansion into Eq. ( |22[ ) and 
comparing the coefficients of powers of s on both sides. 
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The boundary conditions on the probability current must 
also be satisfied by Hk for each k. 

We see that these differential equations may be solved 
iteratively. Fortunately, it turns out that only H 2 {x) 
is needed for our purposes as a short calculation using 
Eq. (19) shows that 


or 


nOO pOQ 

/ dt' mi{t') = — dx dxH^ix) dxi4'{x,\). (26) 

Jo j — OQ 

For potentials that grow unbounded as |x| —)■ oo, the 
probability current must vanish in the limit of large |a:|. 
With these boundary conditions it is not difficult to show 

H2{.x) = ^ 


— OO J a 

{ 9 {x' - Xo) - neq(x', A))^neq(x, A) + 
f dx'/ dx" 

^ —OO ^ a 

(d(x"-xo)-ne,(x",A))l, (27) 


where 9 denotes the Heaviside function and neq(x, A) = 
J^^dx' peq{x',\) is the equilibrium cumulative distri¬ 
bution function. Here, a is an arbitrary real constant. 
Surprisingly, H 2 {x) is independent of a and we will have 
occasion to choose different convenient values for compu¬ 
tational purposes. 

From this result we see that the inverse diffusion tensor 
has the compact form 


Cij(A) = / dv{x,x',x") 


^l3{U(x',X)-U{x,X)-U(x",X)) , 


dxi(l){x,X) dx,(l)[x'',\) 


, (28) 


where we have used the shorthand 
1 


dv{x, x', x") —>■ — 


DZ{X) 


dx 


dx' 


dx". 

O 

(29) 

This expression may be further simplified by observing 
that 

[ dx" aA.0(x",A)e-^^(“"’^) = -Z{X)dx.n,q{x',\) 

J — oo 

(30) 

and 

dxicj^ix, A)e-^^(“’^) = -Z{X)dx.peq{x, X) (31) 

which follow from the definition of the nonequilibrium 
potential cj). These expressions may be used to rewrite 
Eq. (|2^ as 


QjW = - 


ZW 

D 


dx 


f dx' e^^^^'-^^dxmeq{x',X) 

J a 

(32) 


Co(A) = - 


z{\) 

D 


dx 


dxl dxiTIeq(x,X) ) X 


If 


f dx' e^^^^'’^^dx.U,q{x',X) 

J a 

(33) 

lim dxiIleqix,X) [ dx'e^^‘'''''’^'>dxjTleqix',X) = 0, 

x^±oo J ^ 

(34) 


then 


Cii(A) — 


dx 


l9A-neg(x,A) dxiUeqjx.X) 

Peq{x, A) 


, (35) 


which is our desired result. Eq. ( p^ allows us to bypass 
the need for computing correlation functions in order to 
find the inverse diffusion tensor as it is based solely on 
equilibrium information. We will now explore an alter¬ 
nate approach to deriving this formula, and then we will 
apply it to two specific systems of interest. 


III. DERIVATIVE TRUNCATION 
APPROXIMATION 

We relied on linear response theory to arrive at 
Eq. 0- However, the inverse diffusion tensor arises nat¬ 
urally from a first order expansion in temporal derivatives 
of the control parameters as noted in m for harmonic 
potentials. Assuming the probability distribution p{x, t) 
may be well approximated by 

p(x, t) « peq{x, X{t)) + g(x, A(t)) • ^ (36) 

where Q is determined by the Fokker-Planck equation, 
we provide an alternative construction of the inverse dif¬ 
fusion tensor via the so-called “derivative truncation” ar¬ 
gument m- 


Referring to Eq. (10), we assume the nonequilibrium 


probability density has the approximate form 

d\i 

p{x, t) « peq{,x, X{t)) + A(t)), 


(37) 


where Qi(x,X) is to be determined. Substituting this ex¬ 
pression into Eq. (10) and neglecting higher-order deriva¬ 
tives, we see 


dpeq{x,X) 


dX^ 


= D 


' d f d{(3U{x,\j)^ ^ 

& ( — & — »*'*• 

(38) 

Furthermore, since both p(x, t) and Peq{x, A) are normal¬ 
ized probability distributions, we have the constraint 


/ OO 

dx Qi{x, A) = 0. 

-OO 


(39) 
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We may systematically integrate Eq. (381 to obtain a so¬ 
lution that also satisfies Eq. (391 and appropriate bound¬ 
ary conditions. For our purposes, a more expedient way 
of arriving at the solution is to simply state a candidate 
and demonstrate that it satisfies the necessary require¬ 
ments. 

Our candidate is 


and so we must show 


^ -{dxHi[x]\,XQ)dxj<l)(xQ,\))eq,x. (43) 

For unbounded-potential boundary conditions, 
Hi{x-, A, Xo) = d{x — Xo) + neq(x. A). Therefore, 


Gj{x,\) = -{dj;H 2 {x;X,xo)dxj(l)ixo,X))eq,x. (40) 


Here, H 2 is defined by Eqs. (23), (24) and (25) and ap¬ 


propriate boundary conditions as dictated by the physics 
of the problem. The average {■)eq,x applies to the vari¬ 
able Xo and is defined in terms of the stationary state 
probability distribution characterized by A. 

We can quickly establish Eq. (39). Bringing the in¬ 
tegral dx inside of the stationary state average in 
Eq. ( [40| , we see that dx dj;H 2 {x; A, xq) vanishes by 
the fundamental theorem of calculus in the case of an 
unbounded potential at ± 00 . 


Substituting Eq. (40) into Eq. (38), we see that we 
must prove 


= - Dd, [{{^u'{x, A, Xo) 

-f d^H2{x-, A, xo))dx3(t){xo, A))^_^ J. (41) 


(x, A, Xo)r^^ja))(xo, A))eg — (Ja, (()(xo, A))eg _\. 

(44) 

The second term vanishes since peqix,\) is independent 
of Xo and {Oxifpi^o, ^))eq,x — 0. Moreover, 

{^x,XQdxj ^{xq^ X)) eq^X — Peqdxj^ — Peq^ (45) 

establishing the claim. 

We are now in position to relate this derivative trun¬ 
cation approximation to the inverse diffusion tensor ap¬ 
proximation. Recall that 


(>^)a = 



'd\^' 

dt 



(46) 


where 


From Eq. (25) follows 
D [/3f7'(x, \)dxH 2 {x\ A, XQ)+d^H 2 {x-, A, xq)] 


= J_^dx p(x,t)|^(x,A(t)). (47) 

77i(x; A,Xo) 

(42) Using the derivative truncation approximation, 


d\^ f ^ dcf) 

^~~dt J ('^^■^2(2:; A(t),xo) 5 Aj<))(xo,A(t)))g,_x(t)^(a^>A(t)) 

/ oo a / 

dx dj;H 2 {x;\{t),xo)-^{x,X{t)) 

poo poo 

J dt' j dx p{x,t'-,xo)-^{x;X{t)) 

J\j poo 

= -^ cli' {dx-(t>it')dx^(l^iO))eqMt)- 


_dX^ 

dt 

dXi 

dt 


dxiHxoAit)) 

eq.X(t) 

dxj(l){xo,X{t))\ 

/ eq.X{t) 


(48) 


Therefore, we see that the derivative truncation approx¬ 
imation reproduces the inverse diffusion tensor for gen¬ 
eral potentials in the overdamped regime [compare with 

Eq. 0]. 

We note that an expression similar to our Eq. 0 ap¬ 
pears in |53j in which the authors consider anisotropic 
Langevin dynamics in two dimensions. A separation of 
time scales gives rise to slow dynamics for one of the sys¬ 


tem coordinates which can be accurately modeled using 
a Markovian description with a position-dependent fric¬ 
tion. This friction in turn is related to the time integral 
of an autocorrelation function similar to the one used 
here to define the inverse diffusion tensor. A possible di¬ 
rection for future work would be to investigate potential 
connections between the two results. 






















6 



P U{X = 2) 



PU{X= 1.5) 



pU{X=l) 



p U(X = 0) 


rier and waiting a sufficiently long relaxation period [48] . 
After this relaxation period has elapsed the particle has 
equal probability of being located to the left or right of 
the origin. If the particle is found to the left (right) of 
the origin, the memory value is defined to be 1 (0). 

We are primarily interested in optimizing finite-time 
erasure efficiency over cyclic protocols for the classical 
single bit model described above. In m, constraining 
the initial and final probability distributions forced the 
optimal protocols to have jump discontinuities at the end 
points. This was to be expected based on experience with 
optimization in the context of stochastic thermodynam¬ 
ics in general pTHHHl [Ml IM] and erasure efficiency in 
particular [56l El]. These jump discontinuities warrant 
caution when defining thermodynamic quantities such as 
the average dissipated heat ETlESj. 

When classical information is being erased, the differ¬ 
ence in Shannon entropies of the final and initial prob¬ 
ability distributions must satisfy AS = Sf — St < 0, 
which would allow us to define the erasure efficiency 
e = —AS/ {kB{PQ)A) as the ratio of this decrease in 
Shannon entropy to the average heat {Q)\ released into 
the thermal bath (SSJEnj- Taking ks = 1, we see that 


where 


FIG. 1. Continuous erasure protocol. The lefthand well of the 
double-well potential merges with the right and the central 
barrier lowers simultaneously as A decreases from 2 to 0. 


IV. THE ERASURE MODEL 

We consider the following model to represent a single 
classical bit of information: an overdamped Brownian 
colloidal particle diffusing in a one-dimensional double¬ 
well potential in contact with a thermal bath of temper¬ 
ature T [3S| . The wells are initially separated by a 

potential barrier whose height is much larger than the 
energy scale = ksT set by thermal fluctuations, en¬ 
suring stability of memory. Explicitly, we may write the 
potential as 


+ A5 (51) 

is the total average entropy production. Our goal will 
be to minimize the total average entropy over protocols 
while constraining AS. The constraint forces us to op¬ 
timize over protocols with jump discontinuities at the 
endpoints. 

We may express the total average entropy in terms 
of the initial and final probability distributions as well 
as the average work done over the erasure stage of the 
cycle [571IM] : 


= m)A, 


dx p[PU -b Inp] 


T 


(52) 


where r is the duration of the erasure stage. This follows 
from 


U{x,X) = -^log 


ae 


— q;(t —l + A) 


ae 


q;(t—1) 


(l-b (l-be“(^-b)" 


(49) 

where x is a dimensionless spatial coordinate and a ^ 1. 
Initially, A = 2 and there are two distinct wells and a 
central barrier with height governed by a. As A decreases 
to 0, the barrier height diminishes and the left-hand well 
shifts to merge with the right-hand well. 

The system is prepared so that the particle has equal 
probability of being found in either well. This may be 
achieved, for example, by selecting the initial position of 
the particle to be at the midpoint of the potential bar- 


(/3W)A„,3t = f pix,T)[/3U{x,0) - PU{x,t)]. (53) 

Jm. 

Using the approximation Eq. (37l, we see that 

(WA„...=/'<KaAW)(§)Al„(|). (54) 


During the first (erasure) stage, the initial equilibrium 
distribution transitions to a final nonequilibrium distri¬ 
bution in which the system is overwhelmingly more likely 
to have memory value 0. In the second (reset) stage, A is 
brought instantaneously back to its original value while 
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keeping the particle probability distribution constant. No 
heat is generated during this stage. 

Using the calculus of variations, we seek the twice- 
dilferentiable optimizer of dt (^) for the era¬ 

sure stage of the cycle. By considering variations that 
vanish away from the endpoints of the protocol, this op¬ 
timizer satisfies the Euler-Lagrange equation 

dX Ix'^dz ^/az)/r 

vcwm ■ 


where Xi and A/ are determined by the constraints on 
the probability distribution. 

By definition. 


AS = — / dx p{x,T)\n p{x,t) + / da; p(a;, 0) lnp(a;, 0). 

JR JR 

(56) 

Using Eq. (371, we can derive an approximate expression 
for AS, which allows us to compute the endpoints of the 
protocol. The average total entropy generated during the 
optimal protocol is approximately 




(/a? 


+ 0{1/t‘^). (57) 


Since we are neglecting terms of order 1/r^ and higher, 
we need only compute Xi and Xf assuming p{x,t) « 
Peq{x, X(t)) throughout the driving process. For simplic¬ 
ity, we assume that Ai = 2 , A/ = 0. 

From Eq. (35), we have 


C(A) 


1 

TD 


1 - 


tan ^ (^2/[cosh(Q!A) — l] 


^2[cosh(aA) — l] 


(58) 



FIG. 2. Dependence of the leading order term of (As*°*)Aopt 
on a for the erasure model given by Eq. ( |49[ ). Note the satu¬ 
ration at 2 as a becomes very large. 

The dependence of the leading order term of the op¬ 
timal total average entropy on a is displayed in Fig. 
We note that the leading order term grows as a func¬ 
tion of a and appears to approach a limiting value close 
to 2 as a becomes very large. In a similar erasure 


model is considered. There, the potential is piecewise 
constant with central barrier height and “tilt” of the left- 
hand well as control parameters. The leading order term 
of (As*°*)A„pt that model is given in terms of the 
Hellinger distance between the initial and final probabil¬ 
ity distributions: 


-1 2 


Vprir) - Vpr(.0) + 4 x/pi{t) - x/pliO) ■ (59) 


Assuming the particle has overwhelming probability to 
be located in the righthand well after the erasure cycle, 
we see that r(As*°*)A„pt ~ 2.343 for barrier heights of 
physical relevance. (Compare with the current model’s 
saturation value of r(As*°'’)Aopt ~ 2.) In future work, it 
will be interesting to explore whether there might be a 
common cause for the intriguing similarity between the 
saturation values for the leading order terms in the op¬ 
timal average entropy production for these two different 
erasure cycle models. 


V. THE RATCHET MODEL 


We consider an overdamped colloidal particle diffusing 
in one dimension subject to optical trap confinement and 
a tilted ratchet potential m- Mathematically, 

(3U{x, A) = (x - xo)^ - fdFx + pVo tp (x/l). (60) 

Here, I is some characteristic length scale, F is the “tilt” 
of the ratchet, and Vq is the magnitude of the ratchet po¬ 
tential. Furthermore, we choose a single control param¬ 
eter, namely xq, the position of the center of the optical 
trap. 

The model possesses relative mathematical simplicity 
and captures essential physics of chemical processes rele¬ 
vant to the operation of cellular machinery. We view the 
following as the first step towards applying the inverse 
diffusion tensor framework to optimization of nonequilib¬ 
rium processes underlying the functionality of molecular 
motors and other nanoscale biological machines. 

If we consider an isothermal chemical reaction in the 
presence of a catalyst protein {i.e. an enzyme), then 
the reaction can be described by a single reaction coor¬ 
dinate, cycling through a number of chemical states in 
the simplest case m- A suitable working model is then 
an overdamped Brownian particle (reaction coordinate) 
in the presence of thermal fluctuations in a periodic po¬ 
tential. For this reason, we select ip{y) = siii{y). 

If the concentrations of the reactants and products 
are away from their equilibrium ratio, then the catalyst 
molecule will loop through the chemical reaction cycle 
preferably in one direction [50] . In the corresponding 
ratchet model, the periodic potential must be supple¬ 
mented by a constant tilt, i.e. F. 

























Define 


a = j3kP , y = x/l, f = F/(kl) , 
uo{y,yo) = ^{y-yof-fy, e = Vo/{kf) (61) 

so that j3U{x,xo) = ct[uo{y,yo) + e<p(y)] • Assuming the 
strength of the optical trap far exceeds the strength of 
the ratchet potential, e is a small parameter and we may 
apply perturbation theory; i.e. we expand all quantities 
to first order in e and discard higher order terms. 

It is straightforward to show that 


Peqiy,yo) 


Z V 27 r 


1 - ae{<f{y) - {<f}o) 


(62) 

where the subscript “0” indicates an average with e = 0. 


Using Eq. (351, we compute 


Ciyo) ~ (1 - 2ae((^")o). 


(63) 


The equilibrium distribution for the ratchet model 


(Eq. ( |62[ )) has cumulative distribution function 

/ V 

dy'peq{y',yo)- (64) 

-oo 

According to EO], 


e ^ dz < 


w + \/w'^ + - 


(65) 


for -u; > 0. Since ip and all of its derivatives are bounded 
on K, 

lim dygIleq{y,yo) [ dy' dy^U^qiy', yo) = 0. 

y-i-oo 

( 66 ) 

A similar argument shows that 

lim dygIleq{y,yo) [ dy'dy^Hegiy',yo) = 0, 

V^-oo Jp 

. 

and so we may apply Eq. (35) to compute the inverse 
diffusion tensor. 


Using Eqs. (62) and (64), we see that 


[Pva^eq 

Peq 




- (ifiy) - + 2 J" dy' (<p'(j/') - {ip')o) e" t(«|. 

( 68 ) 


From Eq. (35) we have 


72 ( roo py 

C(2/o) « — U + 2ae / dy dy' {(f'iy') - {ip')o) 

L 'J —oo j —OO 


c-'siy'-yo-fY 1^1 — 
a 


(69) 


The integral may be evaluated using integration by parts. 
Therefore, we have 

C(2/o)«^(l-2ae(^")o). (70) 


For ip{y) = sin(j/) (which is mathematically similar to 
the parametric quantron of m), 

C(yo) « ^ (l + 2e e"5);sin(yo + /)) ■ (71) 

Since 

(/31Ue.)A«^"dtC(yo) , (72) 

the minimal dissipation is 


(PWex) Aopt 


(i::SF‘ v^)‘ 


(73) 


We plot the dependence of this expression on / in Fig. 
which clearly illustrates the periodicity in the tilt. This 
periodicity has a natural interpretation. A shift in the 
tilt / by 27r units corresponds to an overall shift in the 
potential energy by a constant amount and a horizontal 
translation upon completing the square. Since constant 
shifts in the potential are physically irrelevant, we would 
expect the optimal average excess work to be insensitive 
to such a shift in the tilt. 


VI. DISCUSSION 


In this paper, we constructed a compact formula 
(Eq. ([nj) for the inverse diffusion tensor for a broad class 
of overdamped dynamics. This is a powerful expression in 
part because the inverse diffusion tensor components are 
expressed entirely in terms of the equilibrium probability 
distribution, and we expect it to increase the number of 
applications for this approach. 

We applied this formula to calculate the optimal dis¬ 
sipation for two model systems. In the first model, we 
considered the erasure of a single classical bit of infor¬ 
mation. The system modeling the storage and erasure 


T 
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FIG. 3. Dependence of minimal dissipation for ratchet model 
on (scaled) tilt / where a = 1, e = 0.1 and i/o(0) = 0 and 
yo(j) = 2. Note that the minimal dissipation is a periodic 
function of /. 


of this bit consisted of a colloidal particle diffusing under 
Brownian dynamics in a double-well potential. In the sec¬ 
ond model, we considered a colloidal Brownian particle 
coupled to an optical trap and a tilted ratchet potential. 
In both cases, the inverse diffusion tensor allowed us to 
predict the optimal time course and produced an explicit 
expression for the minimal dissipation. 

Consistent with previous studies [3711SSH5H], we find 
that the minimal dissipation for each model considered 
here obeys the expected 1/r law. As with the erasure 
model we considered previously the coefficient mul¬ 
tiplying the 1 /r term is fixed by the diffusion tensor, but 
it does not obviously take the same functional form of the 
square of the Hellinger distance between initial and final 
probability distributions. A direct comparison between 
the two erasure models is challenging given that the pre¬ 
vious example was a discrete three state system whereas 
here our system is continuous. Yet, despite these differ¬ 


ences, we found a striking similarity between the satura¬ 
tion values for large, physically relevant barrier heights. 
It is possible that the saturation value for the optimal av¬ 
erage entropy produced during an erasure cycle possesses 
some universality as a function of initial barrier height. 
Future work will address this intriguing possibility. 

Analyzing these erasure models in the regime of low 
to moderate initial barrier height provides another in¬ 
teresting avenue for future development. As described 
in [62], characteristic time scales of dynamical systems 
driven by noise — such as the relaxation time to a steady 
state — follow from the moments of the transition time 
if a specific potential {e.g., Eq. (49)) is assumed. It may 
be possible to apply the framework of [62l |63| to these 
erasure models especially in light of the fact that the 
inverse diffusion tensor may be used to compute the in¬ 
tegral relaxation time for a broad class of systems driven 
by noise [26] . 

We expect further study of these models to generate 
yet more testable predictions for experiments in erasure 
similar to the setup of [64j and the parametric quantron 
based on the Josephson effect in superconductors as 
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